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Abstract 

A more connected world has brought major consequences such as facilitate the spread of 
diseases all over the world to quickly become epidemics, reason why researchers are concentrated 
in modeling the propagation of epidemics and outbreaks in multilayer networks. In this networks 
all nodes interact in different layers with different type of links. However, in many scenarios 
such as in the society, a multiplex network framework is not completely suitable since not all 
individuals participate in all layers. In this paper, we use a partially overlapped multiplex network 
where only a fraction of the individuals are shared by the layers. We develop a mitigation 
strategy for stopping a disease propagation, considering the Susceptible-Infected-Recover model, 
in a system consisted by two layers. We consider a random immunization in one of the layers 
and study the effect of the overlapping fraction in both, the propagation of the disease and the 
immunization strategy. Using branching theory, we study this scenario theoretically and via 
simulations and find a lower epidemic threshold than in the case without strategy. 
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I. INTRODUCTION 


In the last years the complex networks analysis has been focused in no further consid¬ 
ering networks as isolated entities, but characterizing how networks interact with other 
networks and how this interaction affects processes that occurs on top of them. A ^stem 
composed of interconnected networks is called a Network of Networks (NoN) ij-l^. In 
NoN there are connectivity links within each individual network, and external links that 
connect each network to other networks in the system. Very recently physicists have be¬ 
gun to consider a particular class of NoN in which the nodes have multiple types of links 


across different layers called multiplex or multilayer networks 
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Recently, the study of the effect of multiplexity of networks in propagation processes 


such as epidemics has been the focus of many recent researches 
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15|. In Ref Q the 


research concentrated in the propagation of a disease in partially overlapped multilayer 
networks, owing to the fact that individuals are not necessarily present in all the layers of a 
society and this has an impact in the epidemic propagation. For the epidemic model they 
used the susceptible-infected-recovered (SIR) model 17H19| that describes the propagation 


of non recurrent diseases in which infected individuals either die or, after recovery, become 
immune to future infections. In the SIR model each individual of the population can be 
in one of three different states: Susceptible, Infected, or Recovered. Infected individuals 
transmit the disease to its susceptible neighbors with a probability /3 and recover after 
a fixed period of time t^. The spreading process stops when there is only susceptible 
and/or recovered nodes. The dynamic of the epidemic is controlled by the transmissibility 
T = 1 — (1-/5)*’’, which is a measure of the disease virulence, i.e., the effective probability 
that the disease will be transmitted by an infected individual across any given link. At the 
final state of this process, the fraction of recovered individuals R is the order parameter of 
a second order phase transition with a control parameter T. For T < Tc, where is the 
epidemic threshold, there is an epidemic-free phase with only small outbreaks. However, 
for T > Tc, an epidemic phase develops. In isolated networks the epidemic threshold 
is given by Tc = 1/(k — 1), where k is the branching factor that is a measure of the 
heterogeneity of the network. The branching factor is defined as k = (fc^)/ (k), where (A:^) 
and (k) are the second and first moment of the degree distribution, respectively. Since the 
SIR model presents a local tree structure we can employ the branching theory approach 
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In 


lin a generating function formalism 
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2l| that holds in the thermodynamic limit. 


16| the SIR model was studied, with (3 and constant, in a system composed of 


two overlapping layers in which only a fraction q of individuals can act in both layers. 
In their model, the two layers represent contact networks in which only the overlapping 
nodes enable the propagation between layers, and thus the transmissibility T is the same 
in both layers. They found that decreasing the overlap decreases the risk of an epidemic 
compared to the case of full overlap {q = 1). They also found that the critical threshold 
increases as q decreases, and that in the limit of small overlapping fraction, the epidemic 
threshold is dominated by the most heterogeneous layer, this effect could have important 
implications in the implementation of mitigation strategies. 

Motivated by this, in this work we study a disease spreading process in overlapped 
multiplex networks and an immunization strategy for the epidemic spreading. For the 
strategy, we use a random immunization of individuals in one layer of the network. Those 
immunized overlapped individuals will remain immunized in all layers of the network. 


II. EPIDEMIC PROPAGATION PROCESS 


In our model we use an overlapping multiplex network formed by two layers, A and R, 
of the same size iV, where an overlapping fraction q of shared individuals is active in both 
layers. Figure [U^a) shows schematically the partially overlapped network. The dashed 
lines that represent the fraction q of shared individuals should not to be interpreted as 
interacting or interdependent links but as the shared nodes and their counterpart in the 
other layer. 


For the simulation, we construct each layer using the Molloy Reed algorithm 


22 |, 


and we choose randomly a fraction q of nodes in each of the layers that represent the 
same nodes. In our model we assume that the transmissibility is the same in both layers 
because there is only one disease and all individuals in the system spread with the same 
probability. We begin by infecting a randomly chosen individual in layer A. The spreading 
process then follows the SIR dynamics in both layers, the overlapped nodes in both layers 
have the same state because they represent the same individuals. After all infected nodes 
infect their susceptible neighbors with probability (3 in both layers, the time is increased 
in one, and the states of the nodes are updated simultaneously. Note that because there 
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FIG. 1: Partially overlapped multiplex network with layer size iV = 16 and fraction of shared 
nodes q = 0.625. The total size of the network is (2 — q)N = 22 individuals. The dashed lines 
are used as a guide to show the fraction q of shared nodes. Before the spreading dynamics, all 
individuals are in the susceptible stage represented by black circles. 

are shared nodes the branches of infection can cross between the two layers. Thus the 
probability that, following a random link, a node belonging to the infected branch will be 
reached in each layer can be written as, 


n = (1 - q) [1 - Gf (1 - r/x)] + g [1 - Gf (1 - T!a) G„®(1 - T/b)] , ( 1 ) 

/b = (1 - q) |1 - Gf(l - T/b)) + g [1 - Gf (1 - T/b) G „^(1 - T/b)] , ( 2 ) 


where G\^{x) = is the generating function of the degree distribution and 


/02 _ ^max 

'-^1 ~ 2 ^ k = k ^ 


'^min 

„k—l 


Pi{k) kx^ ^ is the generating function of the excess degree distribution in 


layer i, with i = A,B 


2l|. 


Equation ([T]) has two terms, since the probability Ja to expand an infected branch 
following a random chosen link in layer A, can be written as the probability to reach 
one of the (1 — g) non-overlapped individuals and that the branch of infection expands 
through the k — 1 remaining connections of the individual in layer A, combined with 
the probability of reaching one of the q overlapped individuals and that the branch of 
infection expands through the k — 1 remaining connections of the individual in layer A 
and through the k connections of the individual in layer B. An analogous interpretation 
holds for the equation ([2]). 

The solution of the system of equations ([T]) and ([2]) for all T above and at criticality is 
given by the intersection of the curves Ja and fs- At criticality, this intersection can be 
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derived by solving the determinant equation \ J — I\ =0, where I is the identity and J is 
the Jacobian matrix of the system of equations ([T]) and ([2]). The only possibility to have 
a non-epidemic regime is that none of the branches of infection spread, i.e. fA = fB = 0, 
therefore below and at criticality fA = fB = 0. The evaluation of the Jacobian matrix 
•Jij = (dfjdfi)\f, = fr,=n allow us to obtain a quadratic equation for Tc with only one stable 
solution [231 given by, 


_ [{ka - 1) + (kb - 1)] - \/[{ka - 1) - (kb - 1)]^ + Aq^{kA){kB) 

^ 2{KA-l)iKB-l)-2q^kA){kB) ’ 

where k = 1 + l/T^ is the total branching factor of the system and ka^ i^b are the isolated 
branching factors of layer A and B respectively. For g —)■ 0 we recover the isolated network 
result Tc = 1/(ka — 1), which is compatible with our model in which the infection starts 
in layer A and the disease never reaches layer B. In contrast, when g —)■ 1, we hnd 
that Tc = 1/A — i^b)Y + ‘^{kji){kB). Note that Tc(g —)■ 1) < Tc(g —)■ 0). In general, 
Tc decreases with g. This is the case because an increase in the overlapping between 
layers increases the total branching factor, and therefore the total system becomes more 
heterogeneous in degree, i.e., the total branching factor is equal to or bigger than the 
branching factor of the isolated layers. 


III. IMMUNIZATION STRATEGY 


We study a random immunization strategy on the partially overlapped multiplex net¬ 
work. We start by immunizing a random fraction m of individuals in layer A, before the 
epidemic spreading take place. An immunized individual will be immune to the disease in 
all layers, and therefore can not be infected or infect during all the propagation process. 
Note that, due to the presence of the overlapped individuals, in layer B there will be a 
random fraction mq of immunized individuals. 

After immunizing, we spread a disease in the network, starting by infecting a random 
susceptible non-immunized individual in layer A (patient zero). Thus, the probability 
that reaching a node by following a randomly chosen link, it belongs to a branch of 


infection is given by the sw 


distribution in each layer 


,em of equations ([I]) and ([2]), using a node diluted degree 
24| due to the immunization strategy. Thus with the diluted 


degree distribution we have that the branching factor of the diluted layers are. 
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ka = (1 - m) ka 

^ = (1 - qm) Kb , 


( 4 ) 

( 5 ) 


where ka and kb are the branching factor of the original layers respectively. Note that 
the branching factor is rednced dne to the immnnization strategy increasing the epidemic 
threshold and thus hindering the diseases propagation. 


1 

0.8 

0.6 

H 

0.4 

0.2 


Epidemic-Free Phase 


0.01 0.2 0.4 0.6 0.8 1 

q 

FIG. 2: Phase diagram in the plane T — q for the SIR model in the multiplex network, when 
the random immunization strategy is applied, for different values of the immunized fraction m. 
Both, layer A and B, have Erdos Renyi degree distributions with mean values of connectivity 
(/ca) = 6 and (fes) = 4 for layer A and B respectively. Symbols corresponds to the value of Tc 
for different values of m obtained by numerical simulation with layer size N = 10®, while the 
lines denote the theoretical results obtained numerically from Eqs. © and ([2|) using ka and kb 
given by Eqs. o and ©. From top to bottom m = 0.9; 0.7; 0.5; 0.3; 0.1; 0. Above the lines the 
system is in the epidemic phase for each value of m, and below it is in the epidemic-free phase 
where the disease can not propagate. All simulations were done over 10® network realizations. 

In Figure [2] we show the phase diagram in the plane T — g for different values of the 
immunization fraction m. We consider that both layers have Erdos Renyi degree dis- 
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tributions with mean values of connectivity (A;^) = 6 and {ks) = 4 for layer A and B 
respectively, and we use = 1 and kmax = 40 as the minimum and maximum connec¬ 
tivity respectively in each layer. The lines represent for many values of m obtained 
theoretically from Eqs. CQi and ([2]) while symbols denote the numerical simulation results. 
Above Tc there is an epidemic phase and below Tc only outbreaks exists (non-epidemic 
phase). Fig. [2] shows that Tc has different behaviors with q depending on the value of m. 
From Figure [2] we can see a good agreement between the theoretical predictions and the 
numerical simulation results. 

Note that when g = 0 (not shown) the critical threshold corresponds to an isolated 
layer in which the disease starts, i.e. layer A and where the critical threshold is given 
by Tc = — 1). For g —)■ 0 (g = 0.01) the epidemic threshold converges to the 

threshold of the layer with the bigger branchi^ factor, since in this limit the process is 


dominated by the most heterogeneous layer 


16| . We can observe from Fig. [2] that as 


the immunization fraction increases, the epidemic-free phase widens. When m < 0.7 (see 
Fig. |2]) Tc decreases with g owing to the fact that as the overlapping between the layers 
increases the total branching factor of the network increases. However, for m > 0.7 Tc 
increases as g increases. This last effect can be understood taking into account that for 
m > 0.7 layer A is very diluted, thus the disease spreads mostly through layer H, as g 
increases the immunization fraction mg of layer B increases, hindering the propagation 
through that layer. It is expected that for more heterogeneous networks this strategy has 
less impact in the spreading process, due to the fact that the more heterogeneous the 
network is, the more harder it is to dilute with this strategy. 


IV. DISCUSSION 

In this work we study, theoretically and via simulations, an epidemic spreading and 
a random immunization strategy in a partially overlapped multiplex network composed 
by two layer with an overlapping fraction g. We immunize a fraction m of individuals in 
one layer of the network and study how this process affects the propagation of the disease 
through all layers. We found that for g —)• 0 the critical threshold of the epidemic is 
dominated by the threshold of the most heterogeneous layer for all m > 0. We found that 
there is a regime in which Tc decreases with g due to the fact that the total branching 
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factor of the system increases. This behavior stands for m < 0.7, however for bigger 
valnes of m, increases as q increases, hindering the disease propagation. This last 
effect can be understood taking into account that when m > 0.7, layer A is diluted, and 
as q increases the immunization fraction mq of layer B increases, and the effect of the 
immunized individuals in that layer is stronger. 

We can observe from Fig. |2]that as the immunization fraction increases, the epidemic- 
free phase widens. When m < 0.7 we can see that decreases with q owing to the 
fact that as the overlapping between the layers increases the total branching factor of 
the network increases. However, for m > 0.7, increases as q increases, hindering the 
disease propagation. This last effect can be understood taking into account that as q 
increases the immunization fraction mq of layer B increases and for m > 0.7 the effect of 
the immunized individuals in that layer is stronger. Our study suggests that vaccinating 
or isolating only in one layer with the higher propagation capacity, can reduce drastically 
the total branching factor of the network. As a consequence, the epidemic threshold of 
the system increases significantly, reducing the risk of a disease epidemic in the system. 
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